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We discuss a few examples of rich newly developing field of double par- 
ton scattering. We start our presentation from production of two pairs of 
charm quark-antiquark and argue that it is the golden reaction to study 
the double parton scattering effects. In addition to the DPS we consider 
briefly also mechanism of single parton scattering and show that it gives 
much smaller contribution to the cccc final state. Next we discuss a pertur¬ 
bative parton-splitting mechanism which should be included in addition to 
the conventional DPS mechanism. We show that the presence of this mech¬ 
anism unavoidably leads to collision energy and other kinematical variables 
dependence of so-called cr e ff parameter being extracted from different ex¬ 
periments. Next we briefly discuss production of four jets. We concentrate 
on estimation of the contribution of DPS for jets remote in rapidity. Under¬ 
standing of this contribution is very important in the context of searches for 
BFKL effects known under the the name Mueller-Navelet jets. We discuss 
the situation in a more general context. Finally we briefly mention about 
DPS effects in production of W + W~. Outlook closes the presentation. 

PACS numbers: 11.80.La,13.87.Ce,14.65.Dw,14.70.Fm 


1. Introduction 

It is well known that the multi-parton interaction in general and double 
parton scattering processes in particular become more and more important 
at high energies. In the present short review we concentrate on double 
parton processes (DPS) which can be described as perturbative processes, 
i.e. processes where the hard scale is well defined (production of heavy 
objects, or objects with large transverse momenta). In general, the cross 
section for the double-parton scattering grows faster than the corresponding 
(the same final state) cross section for single parton scattering (SPS). 

* Presented at the XXI Cracow EPIPHANY Conference on Future High Energy Col¬ 
liders 
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Fig. 1. SPS and DPS production mechanisms of cccc. 


The double-parton scattering was recognized already in seventies and 
eighties DQ ei si ana mm. The activity stopped when it was realized 
that their contribution at the center-of-mass energies available at those times 
was negligible. Several estimates of the cross section for different processes 
have been presented in recent years [13 [HI [12l [131 (Ill [13 ES ED [18] . 
The theory of the double-parton scattering is quickly developing (see e.g. 

[I3E3EDE21E3E3E3E6J). 

In Ref. [2?J we showed that the production of cccc is a very good place 
to study DPS effects. Here, the quark mass is small enough to assure that 
the cross section for DPS is very large, and large enough that each of the 
scatterings can be treated within pQCD. 

The calculation performed in Ref. m were done in the leading-order 
(LO) collinear approximation. This may not be sufficient when comparing 
the results of the calculation with real experimental data. In the mean¬ 
time the LHCb collaboration presented new interesting data for simulta¬ 
neous production of two charmed mesons [28] . They have observed large 
percentage of the events with two mesons, both containing c quark, with re¬ 
spect to the typical production of the corresponding meson/antimeson pair 
{^DiDj/cTDiDj ~ 10 %). 

In Ref. [23 we discussed that the large effect is a footprint of double 
parton scattering. In this paper each step of the double parton scattering 
was considered in the ^-factorization approach. In Ref. [30] the authors 
estimated DPS contribution based on the experimental inclusive D me¬ 
son spectra measured at the LHC. In their approach fragmentation was 
included only in terms of the branching fractions for the quark-to-meson 
transition c —>• D. In our approach in Ref. [29] we included full kinematics of 
hadronization process. There we showed also first differential distributions 
on the hadron level that can be confronted with recent LHCb experimental 
data [28]. 
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25 years ago Mueller and Navelet predicted strong decorrelation in rel¬ 
ative azimuthal angle m of jets with large rapidity separation due to ex¬ 
change of the BFKL ladder between quarks. The generic picture is presented 
in diagram (a) of Fig. [2j In a bit simplified picture quark/antiquarks are 
emitted forward and backward, whereas gluons emitted along the ladder 
populate rapidity regions in between. Due to diffusion along the ladder the 
correlation between the most forward and the most backward jets is small. 
This was a picture obtained within leading-logarithmic BFKL formalism 
EH 13313311311133 ESI- Calculations of higher-order BFKL effects slightly 
modified this simple picture mmmmmmmmmm leading 

to smaller decorrelation in rapidity. Recently the NLL corrections were cal¬ 
culated both to the Green’s function and to the jet vertices. The effect of 
the NLL correction is large and leads to significant lowering of the cross 
section. So far only averaged values of < cos(nif>jj > over available phase 
space or even their ratios were studied experimentally m • More detailed 
studies are necessary to verify this type of calculations. In particular, the 
approach should reproduce dependence on the rapidity distance between 
the jets emitted in opposite hemispheres. Large-rapidity-distance jets can 
be produced only at high energies where the rapidity span is large. A first 
experimental search for the Mueller-Navelet jets was made by the DO col¬ 
laboration. In their study rapidity distance between jets was limited to 5.5 
units only. In spite of this they have observed a broadening of the 4>jj dis¬ 
tribution with growing rapidity distance between jets. The dijet azimuthal 
correlations were also studied in collinear next-to-leading order approxima¬ 
tion |48| . The LHC opens new possibility to study the decorrelation effect. 
First experimental data measured at yfs = 7 TeV are expected soon jj49 j. 




Fig. 2. A diagramatic representation of the Mueller-Navelet jet production (left 
diagram) and of the double paron scattering mechanism (right diagram). 

The double parton scattering mechanism of W + W~ production was dis¬ 
cussed e.g. in Refs. [Ill 15m I5TI1521155] . The W + W final states constitutes 
a background to Higgs production. It was discussed recently that double- 
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parton scattering could explain a large part of the observed signal [M] . 
We shall also discuss the double parton scattering mechanism of W + W~ 
production in the present paper. 


2. Formalism used in the calculations 


2 . 1 . cccc production 

Let us consider first production of cccc final state within the DPS frame¬ 
work. In a simple probabilistic picture the cross section for double-parton 
scattering can be written as: 

C7 DPS (pp —)• CCCCX ) = <J SPS (pp —> eeX i) • (J SPS (pp —> CcX 2 ). (1) 

2(J eff 

This formula assumes that the two subprocesses are not correlated. At 
low energies one has to include parton momentum conservation i.e. extra 
limitations: x\ + X 3 < 1 and X 2 + 24 < 1 , where x\ and X 3 are longitudinal 
momentum fractions of gluons emitted from one proton and X 2 and X 4 their 
counterpairs for gluons emitted from the second proton. Experimental data 
[55] provide an estimate of a e ff in the denominator of formula ©. In our 
studies presented here we usually take a e ff = 15 mb. 

The simple formula dTJ) can be generalized to address differential dis¬ 
tributions. In leading-order approximation differential distribution can be 
written as 

da 1 da da 

dyidy 2 d 2 pi t dy^dy A d 2 p 2 t 2a e f f dy 1 dy 2 d 2 p u dy 3 dy 4 l d 2 p 2 t 

which by construction reproduces formula for integrated cross section ©• 
This cross section is formally differential in 8 dimensions but can be easily 
reduced to 7 dimensions noting that physics of unpolarized scattering cannot 
depend on azimuthal angle of the pair or on azimuthal angle of one of the 
produced c (c) quark (antiquark). The differential distributions for each 
single scattering step can be written in terms of collinear gluon distributions 
with longitudinal momentum fractions x\, X 2 , x% and X 4 expressed in terms 
of rapidities y±, y 2 , 1 / 3 , 1/4 and transverse momenta of quark (or antiquark) 
for each step (in the LO approximation identical for quark and antiquark). 

A more general formula for the cross section can be written formally in 
terms of double-parton distributions, e.g. F gg , F qq , etc. In the case of heavy 
quark (antiquark) production at high energies: 


da DPS 


n -dggi X 1 ) •^ 2 ) /D) t J, 2)Fgg(.X±X2, Hi, H2) 

2rr 'ff 

da gg ^ c c{x 1 ,x' 1 ,nl)da gg ^ c c(x 2 ,x , 2 , nl) dxidx 2 dx' l dx ' 2 . (3) 







szczurek printed on April 27, 2015 


5 


It is rather inspiring to write the double-parton distributions in the im¬ 
pact parameter space F gg (xi,x 2 ,b) = g(xi)g(x 2 )F(b), where g are usual 
conventional parton distributions and F(b) is an overlap of the matter dis¬ 
tribution in the transverse plane where b is a distance between both glu¬ 
ons in the transverse plane [56]. The effective cross section in ([!]) is then 
l/a e ff = f d 2 bF 2 (b) and in this approximation is energy independent. 

Even if the factorization is valid at some scale, QCD evolution may lead 
to a factorization breaking. Evolution is known only in the case when the 
scale of both scatterings is the same m m m for heavy object, like 
double gauge boson production. 

In Ref. [27] we applied the commonly used in the literature factorized 
model to pp —>• cccc and predicted that at the LHC energies the cross section 
for two cc pair production starts to be of the same size as that for single cc 
production. 

In LO collinear approximation the differential distributions for cc pro¬ 
duction depend e.g. on rapidity of quark, rapidity of antiquark and trans¬ 
verse momentum of one of them m- In the next-to-leading order (NLO) 
collinear approach or in the ^-factorization approach the situation is more 
complicated as there are more kinematical variables needed to describe the 
kinematical situation. In the /ct-factorization approach the differential cross 
section for DPS production of cccc system, assuming factorization of the 
DPS model, can be written as: 

da DPS (pp —> ccccX) 

dy 1 dy 2 d 2 pi ) td 2 p2,tdy3dy4d 2 p 3 jd 2 p 4t t 

1 da SPS (pp —»• ccX i) da SPS (pp —> 00 X 2 ) 

2 a eff dyidy 2 d 2 p u d 2 p2,t dy 3 dy i d 2 p 3 ^d 2 p i)t ' 


Again when integrating over kinematical variables one recovers Eq. ©• 


a eff 


J d 2 b ( T(b )) 2 


where the overlap function 

T(b) = j /(&!)/&-6) dV 


(5) 


( 6 ) 


if the impact-parameter dependent double-parton distributions (dPDFs) are 
written in the following factorized approximation [221157]: 

^i,j(x 1 ,x 2 -,bi,b 2 -, i4, (4) = F i,j(x i,x 2 \nl,i4) f(b 1 ) f(b 2 ). (7) 
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Fig. 3. Production of cc quark and antiquark via fusion of virtual gluons. 


Without loosing generality the impact-parameter distribution can be writ¬ 
ten as 

r(b,xi,x 2 ;f4,(4) = F (xi,f4) F ( x 2 ,lA) F (b;xi ,x 2 ,nl,i4), ( 8 ) 

where b is the parton separation in the impact parameters space. In the for¬ 
mula above the function F(b\ x\, x 2 , p 2 , /t|) contains all information about 
correlations between the two partons (two gluons in our case). The de¬ 
pendence was studied numerically in Ref. mi within Lund Dipole Cascade 
model. The biggest discrepancy was found in the small b region, particu¬ 
larly for large p 2 and/or p 2 . We shall return to the issue when commenting 
our results. In general the effective cross section may depend on many 
kinematical variables: 

<x e ff{xi,x 2 ,x[,x , 2 ,pj,pl) = ( | d 2 b F(b; xi, x 2 , pi, pi) F i b \ x[, x' 2 , pi, pi) 

(9) 

We shall return to these dependences when discussing the role of perturba¬ 
tive parton splitting. 



2.2. Parton splitting 

In Fig. 0]we illustrate a conventional and perturbative DPS mechanisms 
for cccc production. The 2vl mechanism (the second and third diagrams) 
were considered first in [58]. 

In the case of cccc production the cross section for conventional DPS 
can be written as: 


1 1 


1 

167TS 2 


a( 2v2) 


2 &eff,2v2 


dyi dy 2 d 2 p u dy 3 dy 4 d 2 p 2t 


| M(gg ->■ cc)| 2 X 1 X 1 X 2 X 2 
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Fig. 4. The diagrams for DPS production of cccc. 


X D"(xi,x 2 ,l4,l4) D"(x 1 ,x 2 ,l4,l4) 


( 10 ) 


while that for the perturbative parton splitting DPS in a very similar fash- 
sion (see e.g. m) 


a(2vl) = - 


1 1 


dy i dy 2 d 2 pit dy 3 dy 4 d 2 p 2t - 


r| M(gg 


cc)| 2 X\x\x 2 x' 2 


2ff e //,2ui ' rX ‘ 00 167TS 2 

x (d 99 (x , 1 ,x' 2 ,pI,pI)D 99 (x 1 ,X2,pI,iiI) + D 99 (x , 1 ,X2,nl,fj,l)D 99 (x 1 ,x 2 ,li 2 1 ,nlf) 

( 11 ) 


2.3. Four-jet production in DPS 

In the calculations performed in Ref. [59| all partonic cross sections are 
calculated only in leading order. The cross section for dijet production can 
be written then as: 


da(ij —» kl ) 
dyidy 2 d 2 p t 


1 

167r 2 s 2 


x 2fj(x 2 ,p 2 ) I, 


( 12 ) 


where y±, y 2 are rapidities of the two jets and pt is transverse momentum 
of one of them (identical). 

In our calculations we include all leading-order ij —>• kl partonic subpro¬ 
cesses (see e.g. [601 El]). The /t-factor for dijet production is rather small, 
of the order of 1.1 — 1.3 (see e.g. [62, [63]), and can be easily incorporated 
in our calculations. Below we shall show that already the leading-order ap¬ 
proach gives results in reasonable agreement with recent ATLAS [63] and 
CMS [65] data. 

This simplified leading-order approach was used in our first estimate 
of DPS differential cross sections for jets widely separated in rapidity [59] . 
Similarly as for cccc production one can write: 

da DPS (pp —> 4jets A) 


dyidy 2 d 2 pitdy 3 dy 4 d 2 p 2 t 
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Fig. 5. Diagram representing double parton scattering mechanism of production of 
W + W~ pairs. 


y. 

z —'h,ji,ki,h,i2,j2,k2,h 


c dajiiji ->■ fcjZi) da(i 2 j 2 -A k 2 l 2 ) ^ 

a e ff dyidy 2 d 2 p u dy 3 dy^d 2 p 2t 


and partons j,k,l,m = 


whprp c - I 5 if 11 = * 2 J 2 A hh = k 2 l 2 

l 1 if hji 7^ * 2 j 2 V kih ^ k 2 l 2 

g,u,d,s,u,d,s. The combinatorial factors include identity of the two sub¬ 
processes. Each step of the DPS is calculated in the leading-order approach 
(see Eq. (fT2j) L Above y±, y 2 and 2 / 3 , 2/4 are rapidities of partons in first and 
second partonic subprocess, respectively. The pu and p 2 t are respective 
transverse momenta. 

Experimental data from the Tevatron [55| and the LHC E3ES1EZ1 
provide an estimate of a e ff in the denominator of formula (1131) . As in our 
recent paper [ 68 ] we take cr e // = 15 mb. A detailed analysis of a e ff based 
on various experimental data can be found in Refs. [63 ED]. 


2-4■ W + W production 

The diagram representating the double parton scattering process is shown 
in Fig. E] The cross section for double parton scattering is often modelled 
in the factorized anzatz which in our case would mean: 

a W+W~ = e ff a W+ a W- ■ ( 14 ) 

°qq 

In general, the parameter a qq does not need to be the same as for gluon-gluon 
initiated processes < 7 . In the present, rather conservative, calculations we 
take it to be = 15 mb. The latter value is known within about 

10 % from systematics of gluon-gluon initiated processes at the Tevatron 
and LHC. 
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The factorized model (1141) can be generalized to more differential distri¬ 
butions. For example in our case of W + W~ production the cross section 
differential in W boson rapidities can be written as: 

d a W+ S W- _ 1 da^y da w 

dy+dy. a *ff dy+ dy- 

In particular, in leading-order approximation the cross section for quark- 
antiquark annihilation reads: 

= ( x iQi/i( x i^ 2 ) x 2 qj/ 2 ( x 2 ,IJ- 2 ) + x iq i /i(x 1 ,n 2 )x- 2 q j / 2 (x 1 ,fj, 2 )) 

y ij 

x \Mij^ w ±\ 2 (16) 

where the matrix element for quark-antiquark annihilation to W bosons 
(M contains Cabibbo-Kobayashi-Maskawa matrix elements. 

When calculating the cross section for single W boson production in 
leading-order approximation a well known Drell-Yan iv-factor can be in¬ 
cluded. The double-parton scattering would be then multiplied by K 2 . 

3. Results for different processes 

3.1. cccc production 

We start presentation of our results with production of two pairs of cc. 
In Fig. [6] we compare cross sections for the single and double-parton scat¬ 
tering as a function of proton-proton center-of-mass energy. At low energies 
the single-parton scattering dominates. For reference we show the proton- 
proton total cross section as a function of collision energy as parametrized in 
Ref. jUj • At low energy the cc or cccc cross sections are much smaller than 
the total cross section. At higher energies both the contributions danger¬ 
ously approach the expected total cross section. This shows that inclusion of 
unitarity effect and/or saturation of parton distributions may be necessary. 
The effects of saturation in cc production were included e.g. in Ref. m but 
not checked versus experimental data. Presence of double-parton scatter¬ 
ing changes the situation. The double-parton scattering is potentially very 
important ingredient in the context of high energy neutrino production in 
the atmosphere [B1ESIZ2] or of cosmogenic origin [75). At LHC energies 
the cross section for both terms become comparable. This is a completely 
new situation. 

So far we have concentrated on DPS production of cccc and completely 
ignored SPS production of cccc. In Refs. [7(3 :68j we calculated the SPS 
contribution in high-energy approximation m and including all diagrams 
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Fig. 6. Total LO cross section for cc and double-parton scattering production of 
cccc as a function of center-of-mass energy (left panel) and uncertainties due to the 
choice of (factorization, renormalization) scales (right panel). We show in addition 
a parametrization of the total cross section in the left panel. 


in the collinear-factorization approach |68j. In Fig. |7] we show the cross 
section from Ref. [68]. The corresponding cross section at the LHC energies 
is more than two orders of magnitude smaller than that for cc production i.e. 
much smaller than the DPS contribution discussed in the previous figure. 



Fig. 7. Cross section for SPS production of cccc compared to this for standard cc 
production as a function of collision energy. 

In experiment one measures D mesons instead of charm quarks/antiquarks. 
In Fig. [8] we show resulting distributions in rapidity distance between two 
D° mesons (left panel) and corresponding distribution in relative azimuthal 
angle (right panel). The DPS contribuion (dashed line) dominates over the 
single parton scattering one (dash-dotted line). The sum of the two con¬ 
tributions is represented by the solid line. We get a reasonable agreement 
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Fig. 8. Rapidity distance between two D° mesons (left panel) and corresponding 
azimuthal correlations (right panel). 



Fig. 9. Distribution of invariant mass of two D° mesons. 


with the LHCb experimental data [28]. 

Distribution in the invariant mass of two D°D° mesons is shown in 
Fig. (9[ Again a reasonable agreement is obtained. Some strength is missing 
in the interval 10 GeV < M D op,o < 16 GeV. 

At the LHC the cross section for pp —» cc is still bigger than that for 
pp —>• cc.cc [29]. As shown in Fig. [6] the latter cross section is growing fast 
and at high energy it may become even larger than that for single pair 
production. The situation at Future Circular Collider (FCC) is shown in 
Fig. m Now the situation reverses and the cross section for cccc is bigger 
than that for single pair production. We predict rather flat distributions in 
charm quark/antiquark rapidities. The shapes in quark/aniquark transverse 
momenta are almost identical which can be easily understood within the 
formalism presented in the previous section. 
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Fig. 10. Cross section for one c or one c from the cccc final state at the FCC. 



y 


Fig. 11. Rapidity distribution of charm quark/antiquark for yfs = 7 TeV for 
Mi = ™% t , Ml = m 21- 


3.2. Parton splitting 

As described in the Formalism section the splitting contributions are 
calculated in leading order only. In the calculations performed in Ref. [ 58] 
we either assumed fi\ = rn\ t and or p\ = M^ sl and ^ = ^cc 2 - 

The quantity mu is the transverse mass of either parton emerging from 
subprocess i. whilst M c ^ t is the invariant mass of the pair emerging from 
subprocess i. 

In Fig.[lT]we show the rapidity distribution of the charm quark/antiquark 
for different choices of the scale at y/s = 7 TeV. The conventional and split¬ 
ting terms are shown separately. The splitting contribution (lowest curve, 
red online) is smaller, but has almost the same shape as the conventional 
DPS contribution. We can observe asymmetric (in rapidity) shapes for the 
lv2 and 2 vl contributions. 
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Fig. 12. Ratio of the 2vl to 2v2 cross sections as a function of quark/antiquark 
rapidity. 



Fig. 13. Ratio of the 2vl to 2v2 cross sections as a function of collision energy. 


The corresponding ratios of the 2vl-to-2v2 contributions as a function 
of rapidity is shown in Fig. [T2J 

In Fig. [13] we show energy deprendence of the ratio of the 2vl to 2v2 
cross sections. The ratio systematically decreases with the collision energy. 

Finally in Fig. QT] we show the empirical cr e ff, for double charm produc¬ 
tion. Again <x e // rises with the centre-of-mass energy. A sizeable difference 
of results for different choices of scales can be observed. 

3.3. Jets with large rapidity separation 

In Fig. [15] we compare our calculation for inclusive jet production with 
the CMS data [65] • In addition, we show contributions of different par- 
tonic mechanisms. In all rapidity intervals the gluon-gluon and quark-gluon 
(gluon-quark) contributions clearly dominate over the other contributions 
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Fig. 14. Energy and factorization scale dependence of cr e // for cccc production as 
a consequence of existence of the two DPS components. In this calculation cr e ff, 2 v 2 
= 30 mb and cr e ff, 2 vi = 15 mb. 


and in practice it is sufficient to include only these subprocesses in further 
analysis. 







Fig. 15. Our results for inclusive jet production against the CMS experimental 
data m • In addition we show decomposition into different partonic components 
as explained in the figure caption. 

Now we proceed to the jets with large rapidity separation. In Fig. [Tblwe 
show distribution in the rapidity distance between two jets in leading-order 
collinear calculation and between the most distant jets in rapidity in the 
case of four DPS jets. In this calculation we have included cuts for the 
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CMS expriment [39]: y\ , y -2 € (-4.7,4.7), pu,P 2 t € (35 GeV, 60 GeV). For 
comparison we show also results for the BFKL calculation from Ref. [33] . 
For this kinematics the DPS jets give sizeable contribution only at large 
rapidity distance. The NLL BFKL cross section (long-dashed line) is smaller 
than that for the LO collinear approach (short-dashed line). 



Fig. 16. Distribution in rapidity distance between jets (35 GeV < p t < 60 GeV) 
with maximal (the most positive) and minimal (the most negative) rapidities. The 
collinear pQCD result is shown by the short-dashed line and the DPS result by 
the solid line for = 7 TeV (left panel) and -Js = 14 TeV (right panel). For 
comparison we show also results for the BFKL Mueller-Navelet jets in leading- 
logarithm and next-to-leading-order logarithm approaches from Ref. (33]. 

In Fig. [IT] we show rapidity-distance distribution for even smaller lowest 
transverse momenta of the ”jet”. A measurement of such minijets may 
be, however, difficult. Now the DPS contribution may even exceed the 
standard SPS dijet contribution, especially at the nominal LHC energy. 
How to measure such (mini)jets is an open issue. In principle, one could 
measure correlations of semihard (pt ~ 10 GeV) neutral pions with the help 
of so-called zero-degree calorimeters (ZDC) which are installed by all major 
LHC experiments. 


3-4- Production ofW + W pairs 

It was argued that the DPS contribution for inclusive W + W~ could be 
large [53]. Here we partly report results from Ref. [52]. In this analysis we 
have assumed <J e ff = 15 mb as is a phenomenological standard for many 
other, mostly gluon-gluon induced, processes. Similar value was used also in 
other recent analysis [53j where in addition evolution effects of dPDFs were 
discussed. In our opinion the normalization of the cross section may be an 
open issue [52]. Therefore below we wish to compare rather shapes of a few 
distributions. In Fig. [IS] we show two-dimensional distributions in rapidity 
of W + and W~. For reference we show also distributions for 77 and qq 
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Fig. 17. The same as in the previous figure but now for somewhat smaller lower 
cut on minijet transverse momentum. 



Fig. 18. Two-dimensional distributions in rapidity of W + and rapidity of W for 
the DPS mechanism (left), 77 (middle) and qq (right) mechanism for ^fs = 8 TeV. 


components (see a detailed discussion in Ref. [52]). The DPS contribution 
seems broader in the {yw+iVw-) space than the other two contributions. 

In Fig. [T9] we show invariant Mww mass distribution for yfs = 8 TeV. 
The DPS contribution seems to dominate at very large invariant masses. 

How the situation may look at future high-energy experiments at the 
LHC and FCC is shown in Tableland Fig. I2U1 Now the DPS (conservative 
estimation) is relatively larger compared to other contributions. 
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Fig. 19. M ww invariant mass distribution for different mechanism discussed in 
Ref. (52]. 


Table 1. Cross section for W + W production at different collision energies for the 
dominant qq and DPS contributions. 



qq 

DPS 

8000 

0.032575 

0.1775(-03) 

14000 

0.06402 

0.6367(-03) 

100000 

0.53820 

0.03832 


In experiments one can measure charged leptons and not W± bosons. 




Fig. 20. Our predictions for future experiments. 
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Fig. 21. Invariant mass distribution of the W + W~ system (thick solid line) and 
corresponding distribution for the ^ + /i _ system. No branching fractions are in¬ 
cluded. 


Therefore a detailed study of lepton distributions is needed. As en example 
we show (see Fig. ED distribution of invariant mass of charged leptons 
compared with that for gauge bosons. Only a relatively small shift towards 
smaller invariant masses is observed. A more detailed studies are necessary 
to answer whether the W + W~ distribution can be identified experimentally. 
Several background contributions have to be considered. We leave such a 
detailed studies for future. 


4. Conclusions 

We have briefly review some double-parton scattering processes consid¬ 
ered by us recently. 

First we have shown, within a leading-order collinear-factorization, that 
the cross section for cede production grows much faster than the cross section 
for cc making the production of two pairs of cc production very attractive 
in the context of exploring the double-parton scattering processes. 

We have also discussed production of cccc in the double-parton scattering 
in the factorized Ansatz with each step calculated in the /ct-factorization 
approach, i.e. including effectively higher-order QCD corrections. 

The cross section for the same process calculated in the /ct-factorization 
approach turned out to be larger than its counterpart calculated in the LO 
collinear approach. 
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We have calculated also cross sections for the production of DiDj (both 
containing c quark) and DiDj (both containing c antiquark) pairs of mesons. 
The results of the calculation have been compared to recent results of the 
LHCb collaboration. 

The total rates of the meson pair production depend on the unintegrated 
gluon distributions. The best agreement with the LHCb data has been 
obtained for the Kimber-Martin-Ryskin UGDF. This approach, as discussed 
already in the literature, effectively includes higher-order QCD corrections. 

As an example we have shown some differential distributions for D°D° 
pair production. Rather good agreement has been obtained for transverse 
momentum distribution of D° (D°) mesons and D°D° invariant mass dis¬ 
tribution. The distribution in azimuthal angle between both D 0, s suggests 
that some contributions may be still missing. The single parton scatter¬ 
ing contribution, calculated in the high energy approximation, turned out 
to be rather small. In the meantime we checked that 2 4 4 (gg —» cccc ) 
fci-factorization approach leads to similar results as the collinear approach 
discussed here m- 

We have discussed also a new type of mechanism called parton split¬ 
ting in the context of the cccc production. Our calculation showed that 
the parton-splitting contribution gives sizeable contribution and has to be 
included when analysing experimental data. However, it is too early in the 
moment for precise predictions of the corresponding contributions as our 
results strongly depend on the values of not well known parameters a e f f ^ V 2 
and u e ff 2 V l ■ Some examples inspired by a simple geometrical model of 
colliding partons have been shown. A better understanding of the two non- 
perturbative parameters is a future task. 

We have shown that almost all differential distributions for the con¬ 
ventional and the parton-splitting contributions have essentially the same 
shape. This makes their model-independent separation extremely difficult. 
This also shows why the analyses performed so far could describe different 
experimental data sets in terms of the conventional 2v2 contribution alone. 
The sum of the 2vl and 2v2 contributions behaves almost exactly like the 
2 v2 contribution, albeit with a smaller a e ff that depends only weakly on 
energy, scale and momentum fractions. With the perturbative 2vl mech¬ 
anism included, a e jf increases as yfs is increased, and decreases as Q is 
increased. 

We have discussed also how the double-parton scattering effects may 
contribute to large-rapidity-distance dijet correlations. The presented re¬ 
sults were performed in leading-order approximation only i.e. each step of 
DPS was calculated in collinear pQCD leading-order. Already leading-order 
calculation provides quite adequate description of inclusive jet production 
when confronted with recent results obtained by the ATLAS and CMS col- 
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laborations. We have identified the dominant partonic pQCD subprocesses 
relevant for the production of jets with large rapidity distance. 

We have shown distributions in rapidity distance between the most- 
distant jets in rapidity. The results of the dijet SPS mechanism have been 
compared to the DPS mechanism. We have performed calculations rele¬ 
vant for a planned CMS analysis. The contribution of the DPS mechanism 
increases with increasing distance in rapidity between jets. 

We have also shown some recent predictions of the Mueller-Navelet jets 
in the LL and NLL BFKL framework from the literature. For the CMS 
configuration our DPS contribution is smaller than the dijet SPS contribu¬ 
tion even at high rapidity distances and only slightly smaller than that for 
the NLL BFKL calculation known from the literature. The DPS final state 
topology is clearly different than that for the dijet SPS (four versus two 
jets) which may help to disentangle the two mechanisms experimentally. 

We have shown that the relative effect of DPS can be increased by 
lowering the transverse momenta. Alternatively one could study correlations 
of semihard pions distant in rapidity. Correlations of two neutral pions 
could be done, at least in principle, with the help of so-called zero-degree 
calorimeters present at each main detectors at the LHC. 

The DPS effects are interesting not only in the context how they con¬ 
tribute to distribution in rapidity distance but per se. One could make use 
of correlations in jet transverse momenta, jet imbalance and azimuthal cor¬ 
relations to enhance the contribution of DPS. Further detailed Monte Carlo 
studies are required to settle real experimental program of such studies. The 
four-jet final states analyses of distributions in rapidity distance and other 
kinematical observables was performed by us very recently m- 

Finally we have discussed DPS effects in inclusive production of W + W 
pairs. We have shown that the relative contribution of DPS grows with 
collision energy. In experiments one measures rather electrons or muons 
than the gauge bosons. Whether experimental identification of the DPS 
contribution in this case is possible requires a detailed Monte Carlo studies. 
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